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SUMMARY: In this paper we analyze the problem of the total mass and the 
total mass-to-light ratio of the early-type galaxy NGC 4649 (M60). We have used 
two independent techniques: the X-ray methodology which is based on the tem- 
perature of the X-ray halo of NGC 4649 and the tracer mass estimator (TME) 
which uses globular clusters (GCs) observed in this galaxy. We calculated the mass 
in Newtonian and Modified Newtonian Dynamics (MOND) approaches and found 
that interior to 3 effective radii (i? e ) there is no need for large amounts of dark 
matter. Beyond 3R e dark matter starts to play important dynamical role. We also 
discuss possible reasons for the discrepancy between the estimates of the total mass 
based on X-rays and TME in the outer regions of NGC 4649. 

Key words. Gravitation — Galaxies: elliptical and lenticular, cD — Galaxies: kine- 
matics and dynamics — Galaxies: individual: NGC 4649 — Cosmology: dark matter 



1. INTRODUCTION 



The problem of dark matter in early-type 
galaxies is still one of the most important unsolved 
(and, we might add, poorly understood) problems in 
the contemporary extragalactic astronomy and cos- 
mology (e.g. Samurovic 2007, Chap. 1, for a detailed 
introduction). For decades it was almost universally 
considered that the giant ellipticals are paramount 
reservoirs of dark matter, a prejudice partially based 
on the poor understanding of the Faber- Jackson re- 
lation linking the central velocity dispersion and the 
galactic luminosity The general picture which seems 
to emerge with contemporary observational data, 
suggests that interior to ~2-3 effective radii, R e , (in- 
ner regions) dark matter does not play any impor- 



tant dynamical role (e.g. Samurovic and Danziger 
2005). At larger galactocentric distances (beyond 
~ 3 — 4i? e ) it starts to play more important role; 
its contribution, however, varies from case to case 
(e.g. Samurovic and Danziger 2006, Samurovic 2006). 
Overall, there is still considerably smaller amount of 
observational data for elliptical than for spiral galax- 
ies, which makes establishing general results very dif- 
ficult. 

In order to find a typical mass-to-light ratio 
characteristic for early- type galaxies, one can use 
the paper by van der Marel (1991) who found on 
a sample of 37 bright ellipticals that this quantity 
in the B-band is: M/L B = (5.95 ± 0.25)/i 50 ; thus 
M/L B = 8.33 ± 0.35 for h = 0.70 (the value of the 
Hubble constant used in this paper; see Tegmark 
et al. 2004). This sample pertained to inner parts 
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of ellipticals, and therefore we consider the abso- 
lute upper limit for the visible (stellar) component 
M/Lb ~ 9—10. Any inferred mass-to-light ratio 
above ~ 10 in a given region would thus imply the 
existence of unseen (dark) matter there. Note that 
this conclusion is independent on whether the inner 
regions studied by van der Marel (1991) also con- 
tain significant amounts of dark matter (which is 
highly unlikely). We stress that the acceptance of 
large quantities of dark matter in the astrophysical 
community is not unanimous. One possible (and in 
some cases seemingly viable) alternative is MOND 
(MOdificd Newtonian Dynamics), pionereed by Mil- 

grom (1983), which will also be tested in the present 
paper. 

Recently, new kinematical data of NGC 4649 
which extend out to ~ 6i? e , based on the observa- 
tions of globular clusters (GCs), became available 

(Lee et al. 2008). Samurovic and Cirkovic (2008, 
hereafter Paper I) used these data to study the ve- 
locity dispersion profile using the Jeans equation 
and found that the mass-to- light ratio M/Lb ~ 7 
provides a satisfactory fit of the velocity disper- 
sion throughout the whole galaxy which implies 
the scarcity of dark matter in this galaxy. In the 
present paper our aim is to employ other inde- 
pendent techiques to calculate the total mass (and 
mass-to-light ratio) and compare the estimates ob- 
tained using different approaches (both Newtonian 
and MOND). 

The plan of the paper is as follows: In Sec- 
tion 2 we present the basic observational data re- 
lated to NGC 4649 which we use in this paper. In 
Section 3 we calculate the total mass of NGC 4649 in 
Newtonian mass-follows-light approach based on two 
different techiques: X-rays and tracer mass estimator 
(TME). In this section we also apply the virial theo- 
rem in order to calculate the total mass of this galaxy. 
In Section 4 we calculate the mass of NGC 4649 in 
MOND approach (using three different relationships) 
again based on the two afforementioned techiques, X- 
rays and TME. Finally, in Section 5 we discuss our 
findings and present the conclusions. 



2. THE BASIC DATA 



NGC 4649 (M60) is a giant elliptical galaxy 
in the Virgo cluster. It has a nearby companion, 
NGC 4647 (Sc galaxy at 2/5 from the center of 
NGC 4649). The systemic velocity of NGC 4649 is 
^hcl = 1H7 ± 6 kms -1 . Hereafter we use the dis- 
tance based on this systemic velocity - for the recent 
WMAP estimate of the Hubble constant (Tegmark 
et al. 2004; Komatsu et al. 2008), h = 0.70 we ob- 
tain d = 15.96 Mpc; in this case 1 arcsec corresponds 
to 77.5 pc. For the effective radius, as in Paper I, we 
take R e = 90 arcsec (which is equal to 6.97 kpc for 
the distance which we use) , the value taken from the 
paper by Kim et al. (2006). Note that there are dif- 
ferent estimates in the literature, for example, Lee et 
al. (2008) took R e — 110 arcsec and according to the 
RC3 catalog (de Vaucouleurs et al. 1991), R e = 69 
arcsec. Our adopted value is therefore an interme- 
diate one. The discrepancy does not significantly 
impact our conclusions. 

The kinematics of NGC 4649 which includes 
velocity dispersion and symmetric and asymmetric 
departures from the Gaussian is given in Fig. 1 of 
Paper I. The radial distribution of GCs is given in 
Fig. 2 of Paper I: it was shown that beyond ~ 1 ar- 
cmin the surface number density can be fitted with 
a power law: X cx r -7 , where 7 = 1.285 and r is the 
distance to the galaxy center in projection. 

As in Paper I we here use the observations of 
GCs from the study of Lee et al. (2008). This sample 
consists of 121 GCs (83 blue and 38 red GCs) and 
in all our calculations we always deal with the total 
sample in order to have a larger number of GCs per 
each bin. 



3. TOTAL MASS IN 

NEWTONIAN GRAVITY 



In this Section we present the total mass of 
NGC 4649 (and the total mass-to-light ratio in the 
B-band) in Newtonian gravity and in next Section 
we calculate the same quantities in MOND gravity. 

To calculate Newtonian mass we use two dif- 
ferent, independent techniques: X-rays and tracer 
mass estimator (TME). The estimates and compari- 
son are given in Fig. 1 and in Table 1 . In this Section 
we also use scalar virial theorem as a check for the 
estimate of the total mass of NGC 4649. 



3.1 X-RAY ESTIMATE 



The total mass of NGC 4649 interior to the 
radius r based on the X-ray observations is calcu- 
lated using the following equation taken from Kim 
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and Fabbiano (1995): 



M T = 



1.8 x 10 12 (3/3 + a) 
d 



T 



1 keV 



) (lOOO") 



10 Mpc 



Mr. 



(1) 



where the exponent a is related to the temperature 
of the X-ray halo (T ~ r~ Q ) and in our calculations 
taken to be zero, while (3 is the slope used in the 
analytic King approximation model and taken to be 
(3 = 0.50 (Brown and Bregman 2001). The mass-to- 
light ratio in the _B-band calculated as a function of 
radius r is given as: 



1.16 x 10~ 2 10— (3/3 + a) 



T 



1 keV 



(lOOO") 



d 



10 Mpc 



(2) 



where B is the B magnitude of the galaxy inside ra- 
dius r (Kim and Fabbiano 1995). In our estimates 
for the temperature in Fig. 1 we have used the results 
of Randall et al. (2006): T = 0.80 ± 0.05 keV. 

It is important to note that the two relations 
above assume the existence of hydrostatic equilib- 
rium, which is almost universal in interpreting the 
X-ray profiles. Based on the recent work of Dichl 
and Statler (2007) the optical and X-ray isophotes for 
NGC 4649 are very close: for X-rays (based on the 
CHANDRA observations) cllipticity e x = 0.08±0.03 
and position angle P. A.x = 95±27 whereas in the op- 
tical domain e op t = 0.18±0.01 and P.A. opt = 104.0± 
0.4. We interpret this as an indication of hydrostatic 
equilibrium; for NGC 1399 the discrepancy between 
these two quantities in X-ray and optical domain is 
much larger: e x = 0.34 ± 0.04, e opt = 0.10 ± 0.1, 
P.A. X = 179 ± 7 and P.A. opt = 107.4 ± 0.1. The es- 
timate of the total mass of NGC 1399 based on the 
X- rays is given in Samurovic and Danziger (2006). 



3.2 TRACER MASS ESTIMATOR 



Evans et al. (2003) introduced a "tracer mass 
estimator" method which provides an estimate of the 
enclosed mass based on the projected positions and 
line-of-sight velocities of a given tracer population 



(such as GCs in our case) . We assume that the tracer 
population is spherically symmetric and has a num- 
ber density that obeys a power law: 



P(r) = Po 



(3) 



where a is constant, and radius r ranges between r- m 
and r ou t , inner and outer points of the given popu- 
lation, respectively. The parameter 71 = 2.285 (we 
label this parameter 71 to avoid confusion with the 
exponent found in Fig. 2 of Paper I) was determined 
using the surface density of the tracer population be- 
tween r- m and r ou t. Evans et al. give the following 
formula for the mass (supported by random motion) 
enclosed within r; n and r ou t for the isotropic (" iso" in 
formulae below) case ("los" stands for line-of-sight): 



9i 

GN 



ISO \ ^ 2 r> 



(4) 



where Ri is the projected position of the i-th object 
relative to the center of the galaxy, N is the size of 
the tracer population and the constant Cj so is given 
as: 



4 7 i4-a-7i 1- fan/put) 3 ^ 1 

7T 3-7i Mnn/out) 4 ^ 1 



(5) 



This expression is for the case of an isothermal po- 
tential (gravity field is assumed to be scale-free, 
V> = — v^logr, see Evans et al. 2003 for details); we 
have also used the isotropic case (for visible matter; 
for dark matter nothing is assumed) for the sake of 
simplicity, given the small departures from zero of 
the parameter S4 which describes asymmetric depar- 
tures from the Gaussian (see Paper I). 

To obtain the total mass of the galaxy we 
must also take into account the rotational compo- 
nent which is equal to: 



M rot 



G ' 



(6) 



where _R G ut is the outermost tracer projected ra- 
dius in the sample. Therefore, the total mass of 
NGC 4649 is equal to the sum: 



Aftot = M p + M rot . 



(7) 
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Table 1. Newtonian and MOND mass and M/L estimates for NGC 4649 
The table was split in two parts in this preprint. 



r 

(') 
(1) 


M t X ot ay 
10 n M Q 
(2) 


M xr3 -y/L B 

Mq/Lq 

(3) 


M™ E 
10 n M o 
(4) 


M TME /L B 

Mq/Lq 

(5) 


M M ' sim 
M tot 

10 n M o 

(6) 


M M ' sim /L B 

Mq/Lq 

(7) 


<2 


4.3 ±0.6 


4.8 ±0.7 


4.5 ±0.4 


4.8 ±0.4 


3.8 ±0.5 


4.2 ±0.6 


<4 


8.6 ± 1.2 


12.0 ± 1.7 


8.9 ± 1.2 


11.9 ± 1.6 


6.7±0.9 


9.4 ± 1.3 


<6 


12.9 ± 1.8 


20.4 ±2.9 


12.8±2.1 


19.3 ±3.2 


9.0 ± 1.3 


14.3 ±2.0 


<9 


19.3 ±2.8 


32.8 ±4.6 


15.2 ±3.0 


24.6 ±4.9 


11.8 ± 1.7 


20.0 ±2.8 



r 

(') 
(1) 


M M,std 

10 n M Q 
(8) 


M M ' std /L B 

Mq/Lq 

(9) 


M M, toy 

10 n M o 
(10) 


M M ' to y/L B 

Mq/Lq 

(11) 


<2 


4.2 ±0.6 


4.7 ±0.7 


2.9 ±0.4 


3.3 ±0.5 


<4 


7.9 ± 1.1 


11.1 ±1.6 


5.1 ±0.7 


7.1 ±1.0 


<6 


11.0 ± 1.5 


17.4 ±2.4 


6.8 ± 1.0 


10.8 ± 1.5 


<9 


14.2 ±2.0 


24.1 ±3.4 


8.9 ± 1.3 


15.0 ±2.1 



NOTES - Col. (1): radius interior to which a given quantity is calculated, expressed in arc minutes. Col. 
(2): estimate of the cumulative mass based on the X-ray methodology, expressed in units of IO^-Mq for 
T = 0.80±0.05 keV. Col. (3): mass-to-light ratio based on the X-ray methodology, in the £>-band, in solar 
units for T — 0.80 ± 0.05 keV. Col. (4): estimate of the cumulative mass based on the TME methodology, 
expressed in units of IO^A/q for the isotropic case. Col. (5): mass-to-light ratio based on the TME 
methodology, in the B-band, in solar units. Col. (6): estimate of the cumulative mass based on the MOND 
"simple" model (see text for details). Col. (7): mass-to- light ratio based on the MOND "simple" model, 
in the _B-band, in solar units. Col. (8): estimate of the cumulative mass based on the MOND "standard" 
model (see text for details). Col. (9): mass-to-light ratio based on the MOND "standard" model, in the 
_B-band, in solar units. Col. (10): estimate of the cumulative mass based on the MOND "toy" model (sec 
text for details). Col. (11): mass-to-light ratio based on the MOND "toy" model, in the -B-band, in solar 
units. The distance d = 15.96 Mpc was always used. 
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Fig. 1. Total Newtonian and MOND mass (upper panel) and Newtonian and MOND mass-to-light ratio 
(lower panel). In the Newtonian approach we used "tracer mass estimates" and X-ray halo with temperature 
T = 0.80 keV according to the observational data of Randall et al. (2006). In the MOND approach we tested 
three different models: "simple", "standard", and Bckcnstein's "toy" model (see text for details). Additional 
point in the lower panel at 9 arcmin given with the filled triangle corresponds to the MOND estimate of the 
total mass-to-light ratio based on the "toy" model and the TME methodology. The distance d — 15.96 Mpc 
is used everywhere. 



From Fig. 1 it can be seen that the estimates 
based on two physically very different methodologies 
are very similar (especially interior to ~ 6 arcmin), 
thus providing the strong evidence for the increase 
of mass (and the mass-to-light ratio) in NGC 4649 
in Newtonian gravity. The discrepancy between the 



two methodologies and its consequences will be dis- 
cussed below. 
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3.3 VIRIAL ESTIMATE 



4. TOTAL MASS IN MOND GRAVITY 



In order to establish the total mass of 
NGC 4649, we can also use scalar virial theorem for 
a stationary stellar system (Bertin et al. 2002): 



Re 



K 



yc , 



(8) 



where Y* is the stellar mass-to-light ratio in the given 
band, L is the luminosity of the galaxy and a® is the 
central velocity dispersion referred to an aperture ra- 
dius of R e /8. Ky is the so-called "virial coefficient" 
which takes into account the projection effects. Cap- 
pellari ct al. (2006) showed that Ky — 5 ± 0.1 for a 
sample of early- type galaxies at redshift z <~ 0. From 
this equation the total dynamical mass is: 



M. 



dyn 



G 



(9) 



Using the aforementioned value for Ky = 5, ctq — 
230 km s~ , and R e — 6.97 kpc, one obtains the total 
dynamical mass of NGC 4649: M dyn = (4.2 ± 1.3) x 

10 11 M Q . This value takes into account only random 
motions and it does not take into account the rota- 
tional support so it is necessarily lower than the true 
value. Taking into account Eq. (6), we get at 3R e 
M rot (3i? e ) = (0.9 ± 0.7) x 10 n M Q (for v TO t = 141 
km s _1 ). Thus, the total dynamical mass within 3R e 
(= 3'.45) is equal to (5.1 ± 2.0) x 10 n M Q which is 
very close to the values obtained using X-ray and 
TME methodologies at the same galactocentric dis- 
tance (see Table 1). We note that in our estimate 
of the dynamical mass we have used the value of 
do inferred from the kinematics of the GCs. If we 
use the value of the velocity dispersion based on 
the integrated stellar spectra (Fisher et al. 1995), 
(To = 325 km s~ a higher value of the total dynami- 
cal mass (rotational support is assumed) is obtained: 
M dyn = (9.3 ±2.0) x lO n M . This is even closer to 
the values based on the TME and X-ray methodolo- 
gies (see Table 1). 

The role of the rotation in early-type galaxies 
was discussed recently by van der Wei et al. (2008) 
who studied the sample of these galaxies between 
redshift z = 1 and the present epoch z ~ and found 
that for rotating galaxies Ky is possibly <~ 20% 
larger than the canonical value of 5. They used 
Eq. (9) only, in order to establish the total mass of 
the galaxies in their sample. We see that in the case 
of NGC 4649, the contribution of the rotation, equal 
to M rot /M dyn = 0.9 x 10 n M Q /4.2x 10 n M Q w 0.21, 
is close to that given by van der Wei et al., i.e. the 
mass is ~ 21% larger when we take the rotation into 
account, as given above. 



In this Section we calculate the total mass 
of NGC 4649 in MOND gravity using three differ- 
ent formulas: (i) the "simple" MOND formula from 
Famaey and Binncy (2005), (ii) the "standard" for- 
mula (Sanders and McGaugh 2002) and (iii) the 
Bekenstein's "toy" model (Bekenstein 2004). The 
Newtonian acceleration is written as a^v = a/i(a/ao) 
where a is the MOND acceleration and fi(x) is the 
MOND interpolating function where 



Here, a = 1.35 



+0.28 lft 
-0.42 x 1U 



a 
a Q 

-10 



(10) 



m s 2 is a new universal 



constant required by MOND (Famaey et al. 2007). 
The interpolation function fi(a/ao) shows the asymp- 
totic behavior, for a 3> ao, so that one obtains 
the Newtonian relation in the strong field regime, 
and /i — a/ao for a <<C oo- 

As we have shown in Paper I the "external 
field effect" does not play an important role in the 
case of NGC 4649 so we will neglect it in our calcu- 
lations below. 

The MOND dynamical mass, Mm can be 
expressed by using the Newtonian one within a 
given radius r, Mjy, through the following expres- 
sion (e.g. Angus et al. 2008): 



M M (r) = M N (r) x n(x). 



(11) 



The intepolation function can have different 
forms as given below (we refer the reader to Paper 
I regarding our MOND calculations and the details 
related to the variable x and the expressions for the 
circular velocity u c j rc in MOND approach) . We note 
that the Newtonian mass, Mjv is always based on 
Eq. (1) unless stated otherwise. 



A "simple" MOND formula is given as: 
fj,(x) 



1 



A "standard" MOND formula is given as: 
fj,(x) = 



VTTx^ 



(12) 



(13) 



Finally, for the " toy" model of Bekenstein the 
MOND formula is: 



(j,(x) = 



1 + VTT4x~ 



(14) 



The results for the total mass and mass-to- 
light ratio are given in Fig. 1 and Table 1. The un- 
certainties due to the error for the temperature and 
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the errors for the ao parameter are also given in Table 
1. All predictions imply the existence of dark matter 
beyond ~ 3i? e , with the exception of the Bekenstein 
" toy" model based on the TME methodology, which 
is marginally consistent with no dark matter hypoth- 
esis (see Fig. 1). 



5. DISCUSSION AND CONCLUSIONS 



We used the kinematics of 121 GCs and the X- 
ray observations of the early-type galaxy NGC 4649 
in order to examine its mass profile. We per- 
formed our calculations in Newtonian and MOND 
approaches with the goal of establishing the existence 
(or lack) of dark matter in this galaxy. In Paper I we 
found that the velocity dispersion of this galaxy can 
be fitted with M/Lb ~ 7 implying a low amount of 
dark matter, if any. Judging from the approximately 
constant value of the velocity dispersion in the bins 
at various galactrocentric distance from the center 
we found a hint of increasing mass-to-light ratio in 
the outer bins, i.e. we could not exclude the mass-to- 
light ratio as high as 15 in the outermost bin (beyond 
~ 8 arcmin) (see Fig. 3 of Paper I) . 

Our results are as follows: 

1) Interior to ~ 3i? e (=270 arcsec) we find 
no support for the significant amounts of dark mat- 
ter. Both TME and X-ray methodologies predict the 
mass-to-light ratio M/Lb ~ 12 (in the Newtonian 
approach) which is higher than the value obtained 
using the Jeans modelling (see Paper I) but given the 
error bars (see Table 1) we see that the amount of 
dark matter may indeed be very low. In the MOND 
approach the predicted mass-to-light ratio is lower 
(between ~ 7 and ~ 9) which is obviously consis- 
tent with the absence of dark matter interior to this 
radius. The virial estimate of the total dynamical 
mass based on the velocity dispersion inferred from 
the kinematics of the GCs (ao = 230 km s _1 ) is lower 
than the estimate which we obtain using the veloc- 
ity dispersion based on the integrated stellar spectra 
(<7o = 325 km s" 1 ); the latter value is in a very good 
agreement with the estimates obtained using both 
X-ray and TME methodologies. 

2) Beyond ~ 3i? e the situation becomes more 
complex: both X-ray and TME methodologies pre- 
dict rapid increase of the total mass (and the mass- 
to-light ratio) implying the existence of the dark 
component which starts to dominate. Both method- 
ologies provide a very similar prediction: interior 
to 6 arcmin (4i? e ) the mass-to-light ratio becomes 
M/Lb ~ 20 which implies at least 50% of the dark 
matter contribution. At 9 arcmin (6R e ) the predic- 
tions of the two methodologies differ: the X-rays es- 
timate {M/L b "~ 33 is much higher than that based 
on TME, M/Lb ~ 25): note that the estimates are 
nevertheless consistent given the larger error bars 
in the last bin. The predictions of both method- 



ologies imply significant amounts of dark matter in 
the outer regions of NGC 4649 in the Newtonian 
approach. MOND estimates are lower but for the 
estimates based on the X-rays one still needs sig- 
nificant amounts od dark matter in the outer re- 
gions. There is only one case for which we man- 
aged to have rather low mass-to-light ratio in the 
MOND approach: this is a case based on the TME 
methodology and the Bckcnstcin's " toy" formula. At 
9 arcmin M/Lb = 11.3 ±2.2 which is marginally 
consistent with the absence of dark matter (see Fig. 
1). While this marginal consistency may present the 
"last resort" for the MOND proponents, one should 
also point out that all other encouraging results of 
the MOND approach have been obtained using dif- 
ferent models (which we have dubbed "simple" and 
" standard" ) , which in this case require large quanti- 
ties of dark matter, thus obviating the motivation for 
such a radical approach as MOND in the first place. 

3) The question of the discrepancy between 
the X-ray and TME mass estimates remains open. 
One of the possible reasons could be the deviation 
from hydrostatic equilibrium in NGC 4649. We do 
not consider this as a crucial issue in the present case 
(see the discussion at the end of Section 3.1). One 
possible reason can be the lower gas temperature in 
the outer regions of NGC 4649. This could be related 
to the reheating in the central regions, along the lines 
set by the model of evolving cooling flows by Binney 
and Tabor (1995). Note that the last measured point 
in the study of Randal et al. (2006, table 4) is interior 
to <~ 200 arcsec which is well inside 3R ei i.e. within 
the region where all approaches agree on the pre- 
dicted total mass. In order to resolve this question it 
is necessary to have the X-ray observations which ex- 
tend beyond ~ 3i? e and to have more observed GCs 
in the outer bins so as to reduce the uncertainties 
related to the kinematical profile of NGC 4649. 
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